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Abstract We prove that general correlation functions of both ratios and products of characteris- 
tic polynomials of Hermitian random matrices are governed by integrable kernels of three different 
types: a) those constructed from orthogonal polynomials; b) constructed from Cauchy transforms 
of the same orthogonal polynomials and finally c) those constructed from both orthogonal polyno- 
mials and their Cauchy transforms. These kernels are related with the Riemann-Hilbert problem 
for orthogonal polynomials. For the correlation functions we obtain exact expressions in the 
form of determinants of these kernels. Derived representations enable us to study asymptotics 
of correlation functions of characteristic polynomials via Dcift-Zhou steepest-descent/stationary 
phase method for Riemann-Hilbert problems, and in particular to find negative moments of char- 
acteristic polynomials. This reveals the universal parts of the correlation functions and moments 
of characteristic polynomials for arbitrary invariant ensemble of [3 — 2 symmetry class. 



1. Introduction 

Correlation functions of characteristic polynomials for various ensembles of random matrices 
were investigated by a number of authors in a series of recent papers. Keating and Snaith [42[ ], 
Hughes, Keating and O'Connell |35| [36| demonstrated that averages of characteristic polynomials 
over ensembles of random matrices can be useful to make predictions about moments of Ricmann 
zeta function, and other L-functions. These authors consider ensembles of matrices associated 
with compact groups (the simplest case in the family of ensembles of (3 — 2 symmetry class) 
and derive moments of characteristic polynomials J42J. In subsequent papers (Conrey, Farmer, 
Keating, Rubinstein and Snaith Jh], [|) compute more general (autocorrelations or "shifted 
moments") correlation functions of products of characteristic polynomials. 

Brezin and Hikami (and also Mehta and Normand j|4j) considered correlation functions 
of products of characteristic polynomials for an arbitrary unitary invariant ensemble of Hermit- 
ian matrices. This family of ensembles is characterized by the weight exp [— NTtV(H)] in the 
corresponding probability measure ( V(H) is an essentially arbitrary potential function, N is 
the dimension of the matrix H). Using the method of orthogonal polynomials, they found both 
exact and asymptotic (large N) expressions for the correlation functions, and for the positive 
moments. This enabled them to investigate universality of results and find dependence on den- 
sity of states. Namely, the asymptotic expressions were proved to be factorized in product of 
universal and non-universal (ensemble-dependent) parts. It was found that the universal numer- 
ical pre-factors of positive moments of characteristic polynomials coincide asymptotically with 
those for the unitary random matrices obtained by Keating and Snaith jl2| . Thus it was rather 
naturally to expect that these universal pre-factors should appear in the positive moments of 
Riemann zeta function. 
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While Brczin and Hikami compare the positive moments of characteristic polynomials with 
the positive moments of Riemann zeta function, it is similarly worth to compare negative mo- 
ments of characteristic polynomials with the negative moments of zeta function. Note that such 
comparison makes sense only if the degree of universality of the negative moments is established. 
Indeed, it is clear that only universal parts of moments of characteristic polynomials (universal 
pre-factors, for example) may be related with the corresponding moments of zeta function. 

For the negative moments of zeta function a conjecture is available due to the work by Gonek 
|32[ . Fyodorov in j25| performed the calculations for negative moments of characteristic polyno- 
mials for the simplest case of Hcrmitian random matrices with the Gaussian potential function 
V(H) — H 2 /2, known as the Gaussian Unitary ensemble (GUE). The result in |2[| agreed with 
those by Gonek fl32[ . However, a full comparison was still not possible since the universal results 
were unavailable. 

In this paper we both find the negative moments of characteristic polynomials for any unitary 
invariant ensemble of Hermitian matrices and compare them with moments of zeta function. 

A study of the negative moments of characteristic polynomials is further motivated by the 
recent observation by Berry and Keating in [Q. These authors argue that divergences of the 
negative moments could be determined by degeneracies in the spectrum, or clusters of eigenvalues. 
It is an interesting assertion as clusters should be rare events for a random matrix due to the 
level repulsion. Berry and Keating show that the question of whether the influence of clusters 
is dominant is related to that how the negative moments diverge. More precisely, while the 
negative moments of characteristic polynomials are divergent this divergence can be removed 
once we agree that we consider these moments off the real axis (i.e. on the line shifted from the 
real axis by small parameter 5). Then the moments are well defined and, in principle, can be 
computed. According to Berry and Keating scenario it follows that the 2K negative moments 
are proportional to 8~ K as S goes to zero for unitary invariant ensembles of Hermitian matrices. 
In [^|| it was shown that indeed the negative moments diverge as 5~ K for the case of GUE. 
In this paper we prove the universality of this result, i.e. negative moments of characteristic 
polynomials diverge as 8~ K for all unitary invariant ensembles of Hermitian matrices. 

Another important class of correlation functions includes (product of) ratios of characteristic 
polynomials. As is well known those correlation functions can be used to extract more conven- 
tional n-point correlators of the spectral densities (see, for example, ||^, |30) and the references 
therein). As a simple illustration of this statement we compute correlation functions which in- 
clude ratios of characteristic polynomials for an arbitrary unitary invariant ensemble, and then 
reproduce a well known asymptotic result for the two-point correlation function of the resolvents. 

An even more general class of correlation functions are those that include both products 
and ratios of characteristic polynomials (i.e. when the numbers of characteristic polynomials 
in the numerator is different from that in the denominator). These functions provide a very 
detailed information about spectra of random matrices. For this reason such correlation functions 
(together with autocorrelation functions of characteristic polynomials) are pervasively used in the 
field of Quantum Chaos see Jl], ^9) and references therein. For example, they are used for 

extracting generating functions for such physically interesting characteristics as the distributions 
of the "local" density of states and of the "level curvatures" (see e.g. Andreev and Simons [Q 
for more detail). 

For the particular case of the Gaussian Unitary Ensemble (GUE) the large N asymptotics 
of such correlation functions is known. A straightforward way to study the asymptotics is to 
exploit the supersymmetry technique and its modifications [jl], p9fl . For example, in J25| 

one can find asymptotic values for negative moments for the GUE, and Andreev and Simons 
were the first who obtained the asymptotics for other correlation functions. Moreover, Brezin 
and Hikami in their recent paper [|| have made an effort to apply the supersymmetric technique 
to correlations of ratios of characteristic polynomials for more complicated Gaussian Orthogonal 
and Gaussian Symplectic ensembles. 

However, the rigorous application of the supersymmetry technique (and its various modifica- 
tions) is limited to Gaussian ensembles only. Another common disadvantage of these approaches 
is that they are rather robust. While they work well for the computation of asymptotics, in- 
vestigation of the correlation functions for finite size matrices is hardly possible. Moreover, 
those methods hid a nice determinantal structure (revealed in our recent paper f3l| ) of the exact 
expressions and only yielded the asymptotic result in a form of a sum over permutations [jl], p9| . 
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Thus, for understanding the correlation functions of more general (non-Gaussian) ensembles 
a different procedure is required. In the present paper we solve the above mentioned problems 
and provide a unified approach to general correlation functions of characteristic polynomials for 
unitary invariant ensembles of Hermitian matrices. 



2. Statement of the Problem and the Main Results 



Let H denote a N x N random Hermitian matrix which is an element of a unitary invariant 
ensemble (i.e. that of (3 = 2 symmetry class). Introducing TV-dimensional vector x of eigenvalues 
of the matrix H, one defines the ensemble by the eigenvalue density function Pn(x) (see Mchta 
11), 



P N {x) = [Z N ] 1 exp 



A' 



A 2 (£) (2.1) 



The symbol A(x) stands for the Vandermonde determinant, V(x) is a potential function and Zjy 
is a normalization constant. 

A characteristic polynomial which corresponds to the matrix H is defined as Zn [e,H] = 
det(e — H), This object is a building block for constructing various correlation functions of 
interest, such as correlation functions of products and ratios of characteristic polynomials. Our 
main goal is to provide a systematic method for computing all these correlation functions for 
non-Gaussian ensembles, i.e. for the potential V(x) which is more general then x 2 . 

From this end we first consider the correlation functions of the following types: 

F?(k® = (j[ZN[\ j ,H]Z N yM j ,H]J , (2.2) 

^>=(nffeff) H , m 

where the averages are understood as integrals with respect to the measure djj,(x) = Pj^{x)dx 1 



with P/v(x) being defined by the equation (2.1). If the components of vectors e and u> have 



nonzero imaginary parts the correlation functions above are well defined. 

Our first result is that each of the above correlation functions are essentially governed by two- 
point kernels constructed from monic orthogonal polynomials iTk(x) and their Cauchy transforms, 
ftfe(e). The monic polynomials, itj(x) = x J + . . ., orthogonal with respect to the measure dfi(x) — 
e~ NV ^ x 'dx, are defined by 

n k {x)Tr m (x)e~ NV( - x ' > dx = c k c m S km (2.5) 



and their Cauchy transforms are determined in accordance with the following expression 
Zm J x — e 

The correspondence between the types of correlation functions and different kernels obtained in 
this paper is summarized in Table 1. 

One of those kernels, the kernel Wi^n+k{\ A*)> i s weu known in the theory of random matrices. 
It is related to the familiar kernel K^(X, fi) which is known to determine completely the n-point 
correlation functions of eigenvalue densities as well as spacing distributions between eigenvalues. 
The kernel Kn{\^) is defined by 

K N (x,y) - 7jV ~ 1 e -f n*) Mf>[Wj/) -KN-i(x)*N{y) c -Kv( v ) (2 ?) 

' 2ni x — y 
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Correlation function 


Kernel 


•^F(A,A) 


TTA / \ .A 7riv + A-(A)7rj V+K _i(/i)-7rjv+ji:-i(A)7rjv+K(M) 




WYr iv fc u) ^ h N (e)ir N - 1 (ii)-h N - 1 {6)Tr N (jj,) 






Wiii,n-k(£,u) — e _ u 



Table 1 . Correlation functions and kernels 



where 

27TZ . „. 

7«-i = -~s— (2-8) 

C n-1 

However two other kernels Wii t w(e, fi), Wih^n-k^-, 1 ^) responsible for the correlation functions 
of characteristic polynomials have not been previously considered, to the best of our knowledge. 

Apart from the correlation functions discussed above we investigate also the correlation func- 
tions containing non-equal number of characteristic polynomials in the numerator and the de- 
nominator, such as 

lY[Z N \m,H\\ 

^ M (e,A) = (^ ) , 0<M<K (2.9) 



iv ^>p; \ m 

\Y[Z N [e 3: H]l 

1 H 



and 




^ M (e,A) = (^ ) , 0<K<M (2.10) 

\J[Z N [e h H]l 

H 

where K + M is an even number, i.e. K + M = 2L. We reveal that these functions can be 
expressed in terms of determinants of size L x L. The entries of the determinant for the func- 
tion Tfy M (e, p,) are kernels Wu^n-m+l and Wi^n-m+l while the entries of the corresponding 
determinant for the function Ty' M (e, ji) are the kernels Wjj^-m+l an d Wjjx,n-m+L- 

It is important to note that all the kernels introduced above belong to the family of the so- 
called "integrable" kernels, i.e. they correspond to integrable operators first distinguished as a 
class by Its, Izergin , Korepin and Slavnov |3£| |3£| [57j. According to that theory the kernels 
corresponding to the integrable operators are defined as follows: Let E be oriented contour in C. 
An operator L acting in L 2 (T,, C) is called integrable if its kernel has the form 

M 

L(z,z') = ^ ; (2.11) 

z — z 

for some functions fi,gj] i, j = 1, ■ ■ ■ , M. The formalism of integrable operators is described by 
Deift in Q. 

Our kernels obviously satisfy the above definition. This is an important observation as in- 
tegrable kernels have a variety of useful properties that can be exploited (see a review paper 
by Deift, Zhou and Its [^3[). For example, it is known that the Fredholm determinants of inte- 
grable kernels satisfy non-linear differential equations. In particular, Tracy and Widom [pl||-[p4| 
obtained differential equations (of the Painleve type) for the Fredholm determinants associated 
with the kernel K]y(x,y). In the present paper we use another important property of inte- 
grable kernels, namely, their relation with Riemann-Hilbert problems. In a similar way as the 
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Riemann-Hilbert technique was applied to the kernel K^{x,y) to demonstrate its universality 
in the Dyson's scaling limit (see fl^j , |l3f and ^)) we exploit the Riemann-Hilbert approach 
to find Dyson's scaling limit of kernel functions Wi^n+k^, h), Wn^i^, ^) an< ^ Wiii,n—k ( £ j UJ )- 
As a result we obtain associated universal kernels summarized in Table 2. 



Finite N kernel functions 


Associated limiting kernels 






Wjj,iv(e, fi), 




fs!^, 3mC>0 


3m e / 








S/w(C -»?) = < 


^z - , 3m C > 0, 3m 7/ < 
3m C < 0, 3m 77 > 
k otherwise 


3m e ^ 0, 3m u> ^ 







Table 2. Finite and associated limiting kernels 



The representation of the correlation functions in terms of determinants of the kernels (see section 
4) enables us to give explicit asymptotic formulae for all five correlation functions of characteristic 
polynomials discussed in the text above. We give a summary of these results below. 




Define A"-dimensional vectors, x = (x, . . . , x), where x belongs to the support of the equilibrium 
measure for the potential function V(x) (see section 5 for the definitions), ( = (Ci, ■ • ■ , (k) and 
fj = (771, . . . , 1)k)- Then for the correlation function of products of characteristic polynomials we 
obtain 

Tf{x + C/Np(x),x + f)/Np(x)) 



e KNV {x)lN p(x) f 2 e a{x) V il+m) 

= [cn? K [ Z A " det [MC, ~ (2-12) 

A(C)A(7i) 

where p(x) stands for the density of states, and we have introduced the notation 

a{x) = ^f> (2.13) 

K(£ f,\ - / TT z nIvj> h ] 



2.2 Dyson's Limit for Tfj(l, A) = \ TI 



Z N [e 3 ,H] 



H 



Let x, C, 7} be K -dimensional vectors. Assume that the components of £ have non-zero imaginary 
parts. Then we find 

Tfj(x + C/N P (x),x + rj/Np(x)) (2.14) 
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K(K-l) -a(x)E(Cl-»7l) A(C,7?) nwr <C /A M 

(-) 2 e 1 A o,iw det (FMCi ~ *?j)]l<iJ<JC 



A2(C)A2(f}) 

2.S Dyson's limit of T? u (w,Cj) = (n ^r^-g^^g; 



if 



It is convenient to introduce 2i-T dimensional vector, e = (w, Co). The new coordinates appropriate 
for investigation of Dyson's asymptotic limit of the correlation function are defined so that 
e = x + Q/N p{x), dim.x = dim^ = 2K. Here the vector x has 2K equal components x, and x 
belongs to the support of the equilibrium measure for the potential function V(x). As for the 
components of the vector (, we assume that they have non-zero imaginary parts. We find 



Tf u {x + (/Np(x)) = (-) K [ lN \ K [Np(x)] K2 e - KNV ^\ 



-a(x)Y.«l+Cl + K) 



I det [8 HI - CkU+k))] k , j<k 



(2K)\ 



7T G S 2 



A(Ctt(1)! • ■ • i Ctt(K))^(Ctt(K+1), ■ • ■ , Cr(2if)) 



(2.15) 



24 Dyson's limit for Tf v (e, ft) = ( ± ) , K > M 



IlZ N [e j ,H] i 



H 



In this case new coordinates are introduced so that e = x + (/Np(x), fi = x + fj/Np(x). It is 
clear that dim £ = M, dim fj = K, 5m ( ^ 0. We define 2L = K + M (i.e. we consider the 
correlation function of an even number of characteristic polynomials) . With these definitions we 
find 

Tf v (x + C/Np(x),x + fi/Np(x)) 



( _)^^ [ CN ] K - M [Np{x)] {L ~ M ^ e ^-MM*) 



(2.16) 



r m k 

-a(x)\Yj Ci-E m 
x e L 1 1 



A(Ci, ■ ■ • ,(m;Vl-m+i, ■ ■ -,Vk) 



x det 



A 2 (C)A 2 (r7 L _ M+ i, _ _ _ iVK )A{rn, . . . ,r) L - M ) 
§//(& - Vl-m+i) ■ ■ ■ S//(Ci - Vk) 



§// (Cm — Vl-m+i) 
$i(m - Vl-m+i) 



§//(Cm - 



S/(r/z,_M - Vl-m+i) ■ ■ ■ §i{vl-m - Vk) 



2.5 Dyson's limit for Ty (e, p) = ( -jj 

\UZN[e 3 ,H], 



H 



K <M 



We introduce new coordinates x, (, fj as in the previous case. Let 2L = K + M. Then the 
Dyson's limit of the correlation function Ty (x + (/Np(x), x + fj/Np(x)) is 



^(x + C/Np(x),x + fj/Np(x)) 



M(M-l) ,JK-I r., , 

(-) ^ [ 7Ar ] [Np(, 



i(L—K) 2 N{L-K)V(x)_ 



1 1 

A 2(7)) M! 



(2.17) 
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x dct 



x ( L i i J > 



A(?7l, . . . ,T] K -X^l + M-Ky . . . ,Ctt(M)) 



n ft M A2 (C(i+ ■ ■ ■ , Cr(M))A(Cr(l), ■ ■ ■ , C(M^)) 



(C ( i + «_^) - m) 



§77/(^(1) - Cr(l+ flf - K ))' 



S///(C(M^) - C(i + JM_J1; 

2.6 TTie average of the resolvent 



§Il(Cn{M) - Vi) 
Sz/(Ctt(M) - ??a) 

§///(C7r(l) — Ctt(M) 
§JJj(C 7r ( M - g ) _ Ctt(M)) 



An interesting observation is that the non-universal functions emerging in above expressions for 
the correlations of characteristic polynomials can be expressed in terms of the large N limit of 
the averaged resolvent defined as 

R+(x) = (2.18) 



x — H 



H 



Here we assume that the parameter x has an infinitesimal positive imaginary part. It is then 
straightforward to observe that R^ix) is expressed in terms of the correlation function Fjj (e : p). 
We have therefore the following expression for the large N limit of the averaged resolvent: 

NV'(x) 



R+(x) = inNp(x) 
attributing a particular meaning to the non-universal factors: 

, , 1 i- i _ . -i . . 5Re \Rt(x)] 

Np{x) = -3m JJ+ x) , a(x) = -tt - r +\ A 

2.7 Universality of ^^(e, p) at the center of the spectrum 



(2.19) 



(2.20) 



Equation ( [2.14 ) implies that the following theorem is valid 



Theorem 2.1. Assume that a{x) = V'(x)/2p(x) = (the center of the spectrum, for example), 
where V(x) is the potential function, p(x) is the density of states and x belongs to the bulk 
of the spectrum. Then correlation functions of ratios of characteristic polynomials of random 
Hermitian matrices are universal in the Dyson scaling limit. 



Remark. It can be observed that for the Gaussian case V(x) — x 2 /2 our results for the correlation 
functions (at the centre of the spectrum x = 0) are reduced to the formula by Andreev and Simons 
[0 . A detailed derivation of the large N asymptotics for the Gaussian case can be found in . 

2.8 Negative moments of characteristic polynomials 

Positive moments of characteristic polynomials are determined by the following asymptotic ex- 
pression (Brezin and Hikami |||): 

(Z% K [x,H]) H = [c N f K e KNV ^ [Np(x)f 2 T+ (2.21) 

where 

A"-l 

t a- = n i -/( i + K )- ( 2 - 22 ) 
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T J- is a universal coefficient which also appears in positive moments of zeta-function. 

In the present paper we, in particular, obtain an asymptotic result for the negative moments 
of characteristic polynomials Ad^ N (8). These moments are defined as 



M« tN {8) = {Z N K (x + ,H)Z N K (x-,H)) H 



where 



r± — 



x± 



1 6 



x and 6 are real parameters and p(x) is the density of states. 



2Np{x) ' 

large TV limit we find that the negative moments behave asymptotically as 



(2.23) 
For the 



[2wf [c N 



-2K e ~KNV(x) 



Np{x) 



A 



Formula (2.24) should be compared with Gonek's conjecture 
Riemann zeta function, which states that 




c 



logT 



it 



-2K 



dt 



logT 



A 



(2.24) 

for the negative moments of 
(2.25) 



A similarity between (2.24) and ( 2.25| ) becomes apparent if we put 
Np{x) = logT 



(2.26) 



in accord with the known expression for the mean density of Riemann zeroes. The pre-factor 
[27r] [cn] 2K in front of the exponent in equation ( 2.24 ) is not universal, and as such it is 
irrelevant for the comparison with the moments of zeta function (e.g. for the Gaussian case it is 
equal to e NK in the large N limit) . On the other hand, the analogue of the universal coefficient 
Ti for the negative moments is given by 

T K = 1 (2.27) 

as is immediately evident from the formula ( [2.24 ). In other words, the universal coefficient for 
negative moments of characteristic polynomials which should also appear in the negative moments 
of zeta function is 1. 

Furthermore, we can see from the expression (2.24) that the negative moments diverge at 
5 — > as 8~ V ( K \ with the exponent v(K) being equal to K 2 . This fact fully agrees with the 
behaviour conjectured by Berry and Keating M for all unitary invariant ensembles of (3 = 2 
symmetry class (see the discussion in the Introduction). 



3. Lagrange Interpolation Formula and Identities for Characteristic Polynomials 



In this section we discuss some consequences of the Lagrange interpolation formula (see Szego 



50 [). The o btai ne d rel ations enable us to derive exact expressions for the correlation functions 



M) 

(U)-S S (HI) 

Let x±,X2, ■ ■ ■ ,xm be eigenvalues of the matrix H. Let us associate with the characteristic 
polynomial Zn[c, H] of the matrix H 



Zn[c,H] = (e- Xi)(e - x 2 ) ■ ■ ■ (e - x N ) 
the fundamental polynomials of the Lagrange interpolation: 

Z N [e,H] 



Z' N [x v , H](e - x v ) 



v = l,2, 



From equation ( |3.l| ) It is easy to observe that 
Z' N [x„,H] = \\_(x u - Xj), v = l,2,.. 



,N 



(3.1) 
(3.2) 

(3.3) 



In particular, the equations ( p. 2] ) and (3.3) imply that the fundamental polynomials of the 
Lagrange interpolation have the followin property: 



(3.4) 
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As each polynomial P(x) of degree N — 1 is determined uniquely by its value in N points, we 
have 

P(x) = P( Xl )h(x) + P(x 2 )l 2 (x) + ... + P(x N )l N (x) (3.5) 
From the expression ( |3.5[ ) it follows that 

N 

v=\ 
N 

e = yixulvje) (3.6) 
v=\ 

JV 



We immediately conclude from the above expressions and the equation ((3[|) that the following 
algebraic identities must hold: 

N 

f £ f x*/Z , N [x v ,ff\ = 0, 0<K<N-2 (3.7) 
i/=i 

and 

f K JL T K 1 

7T1 m = J2 rf\> VK = 0,...,N-1 (3.8) 

Z N [e, H\ ^-J e - x v Z' N [x u , H\ 

With these equations in mind it is not difficult to obtain a representation for the Cauchy transform 
hjsr-x(e) of monic orthogonal polynomial 7T/v-i(x) defined by formula ( |2.6| ) in terms of a multi- 
variable integral (for a derivation see [j3l|): 



1 f TT ~ N ? y(:Ej) 

= 75—/ (e-^) e 3-1 A (11, . . . ,x N )dxi . . .x N (3.9) 
The right hand side of this formula can be looked at as the average of 2^ 1 (e, H) taken over the 



ensemble of unitary invariant Hermitian matrices. In other words, equation (3.9) implies that 

(Z^ 1 [e,H}) H =j N _ 1 h N _ 1 (e) (3.10) 

In what follows we compute more complicated correlation functions of characteristic polynomials. 
Those correlation functions include products of characteristic polynomials both in the numerator 
and the denominator. The algebraic identity which enables us to average the product of the 
characteristic polynomials in the denominator is 

M JV-M 
w 



n 



\Z N [ei,H] 



E 



tt x a(i) M j A(x CT(1) , . . ■ ,x tr(M ))A(a: ff ( M+ i), . . . ,x CT( jv)) (3 11) 



where cr € Sn/Sn-m X Sa/, Sat is the permutation group of the full index set 1, . . . , AT, whereas 
Sm is the permutation group of the first M indices and Sjv-a/ is the permutation group of the 



remaining N — M indices. Identity (3.11) was proved in j3T| and follows as a consequence of the 
Cauchy-Littlewood formula jmj 



M N 



YiYv- 1 ~ Xiy j) 1 = z2 s x{ x i^--^N)sx(yi,---,yM) (3.12) 
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and the Jacobi-Trudi identity EtJ : 



det ( x 



s x (x 1 , . . .,x N ) 



(3.13) 



A(xi, . . .,x N ) 

where the Schur polynomial s\(xi, . . . ,xn) corresponds to a partition A, and the indices i,j take 
their values from 1 to N. 



4. Finite Correlation Functions 



A" 



4-1 Correlation function Tj (A, A) = ( II %n [Aj, H] Zn [fij, H] 

\ i / H 

The correlation function of products of characteristic polynomials Tf- (A, A) was investigated in 
detail by Brezin and Hikami H, Q| . For finite size N these authors have demonstrated that the 
correlation function could be rewritten in a determinant form. Namely, 



^f(A,A) = 



1 



A(A,A) 



det 



TTjv(Ai) 7Tat + i(Ai) 

7Tjv(A 2 ) 7T7V+l(A 2 ) 

7rjv(Aic) ir N+ i(\ K ) 

n N (fll) 7T A r +1 (^ 1 ) 

7TAf(^2) 7TAf+l(^2) 



7I"AT+2A'-l(Ai) 
7T7V+2A-l(A2) 

7rw+2A-i(A^) 

KN+2K-l(fl>2) 
7I"A+2A-i(^a) 



(4.1) 



The same correlation function Tj (A, A) has also an alternative representation in terms of a 
determinant of a kernel constructed from monic orthogonal polynomials. 



Proposition 4.1. |7jj: The correlation function of products of characteristic polynomials is 
governed by a two-point kernel function constructed from monic orthogonal polynomials, 

det[W I ^ N+K (X llH )] 1 ^ K (4.2) 



Tf (A, A) = 

A(A)A(A) 

where the kernel Wi t N+K{x,y) is given by the formula 

A — fi 

The constant Cn,k can be expressed in terms of the coefficients 7; defined by equation $2.{\) 



Cn,k — [cn+k-i] 



N+K-l . ,R 

-2K TT („ \2 _ [lN+K-l\ 



n ^ 



YfN+K-1 

11a 



(4.4) 



11 



Proof. To prove formula (4.2) we observe that the correlation function J-^(X,jd) can be repre- 
sented as the integral 



J N 



d N xe 



-N Y, V(xi) 



A(A,x)A(A,x) 



(4.5) 



A (A) A (A) 

This integral can be evaluated using the method of orthogonal polynomials. Namely, we rewrite 
the Vandermonde determinants as determinants of monic orthogonal polynomials. Then the 
product of the Vandermonde determinants in the integrand above can be rewritten as a sum 
over permutations, i.e. 

A(X,x)A(p,,x) 



ffjp £ Sn + k 



R-H APPROACH FOR CORREL. OF CHARACTERISTIC POLYNOMIALS 



11 



(4.6) 



x 7r p(i)-i(A 1 i) • • ■ 7r p(K)-i(A*if) 



x K p{K+i)-x{xi) ■ ■ ■ k p (n+k)-i{xn)] 

We insert the above formula into the integrand of !Ff(X,jl) and integrate over the variables 
Xi,...,xjf. The orthogonality of monic polynomials leads to the expression 

~N+K-1 



^f(A,A) = 



7 -l 

N 



n 



A(A)A(A) 



(7,P £ Sn + K 



?p(i)-i(Mi) • • ■ 1 p (k)-i{^k) 



(4.7) 



X ^CT(_ff+l)p(i<'+l) ■ ■ • ^(K + ATJpfK + AT)] 

Here we introduced the polynomials qi(x) = c^ 1 ni(x) normalized with respect to the measure 
d(i(x) = exp (— NV(x)). The sum in the equation above can be further transformed to a deter- 
minant (see Appendix A) and we end up with the following expression 

~N+K-1 



Fi (A, A) 



n c 

o 



Z N A(A)A(A) 



det 



N+K-l 

E 





9i(Ai)<7;(Mi) 



(4.8) 



l<i,j<K 



Applying to (4.8) the Christoffcl-Darboux formula (see, for example, Szego |30j) we recover the 
expression (4.2). □ 



K 



4-2 Correlation function !Fjj(e, (i) — (f\ ^^gj 

Here we derive an exact formula representing Tf\(j±, e) as a determinant of the kernel Wu^ie, m)- 



Proposition 4.2. Let 3m tj ^ 0, j = l,...,K. Then the correlation function of ratios of 
characteristic polynomials is determined by a two-point kernel constructed from monic orthogonal 
polynomials and their Cauchy transforms. More precisely, the following formula holds 



k A(e,/1) 
A 2 (e)A 2 (/i) 



det [W IItN (e i ,iJ j )} 1 



<i,j<K 



where the kernel Wji^(e, [i) is given by 



e-fj, 



(4.9) 



(4.10) 



and the constant jn-i is defined by the equation 



Proof. We propose a "reduction procedure" . The idea is to reduce computation of the correlation 
functions containing ratios of characteristic polynomials to that of the correlation function which 



contains only products of characteristic polynomials . Namely, we exploit the identity (3.11) to ex- 
press the denominator, Yli ^i^^h as a sum over permutations. T fUe, fi) is a multi- variable 
integral with the measure defined by the eigenvalue density function ( [2.l| ), so the integrand is 
symmetric under permutations of the variables of the integration. (Recall that x\, . . . ,Xn de- 
note eigenvalues of the Hermitian matrix H , dim H = N). It means that each permutation gives 
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the same contribution to the correlation function. The total number of those permutations is 
(n-^k)\k\ • ^ e then find that 



N\ 



(N - K)\K\ 



K 

n 



K-N 



(4.11) 



^ - — T* ■ 



i,3'=l 



A AT 

n n 

,i=l s=A+l 



The next step is to decompose the integration measure in accordance with the following expression 
for the eigenvalue density function 

" K N 

n n (*i-*.) 2 



P (N \xi, . . . ,xn) 



Z 



N-K^K 



1=1 s=K+l 



xPW( Xl ,...,x K )p( N - K Xx K+1 ,...,x N ) 
which allows one to rewrite the correlation function as 



(4.12) 



AT! 



'N—K 



K 

n 



K-N 



{N-K)\K\ Z N 
J dx\ . . . cLxk A 2 (xi, . . . , xk 



(4.13) 



K 

n 

«j=i 



e -NV{ Xi ) x N-K {H _ Xi) 



Y[Zn-k[^i,H}Z n ^ k [xi,H] ) 
\ i In 



Here dim H — N — K, W = H. We then observe that the original integration over TV variables 
is replaced by an integration over K variables. Moreover, we notice that the correlation function 



of products of characteristic polynomials emerges in the integrand of the formula (4.13). Then 
the equation ( |4.2j ) yields 

(f[Z N _ K [ R ,H}Z N _ K [x h H}^ = (4.14) 



Cn- 



K,K 



which leads to essential simplifications in the expression for jj): 

N\ Z N _ K 



— K,K 



{N-K)\K\ Z 



N 



K 



K-N 



/A(A) 



J dx\ . . . d X K A(a;i, . . . , Xk ) 



(4.15) 



det [W / / i Ar(^ i ,X :( )] 1 



<i,j<K 



The first two terms in the integrand can be further rewritten as K x K determinant, 



A(x u 



,x K ) 



K 



n 

2J = 1 



= (-) i ^det [/ 1 (^)] 1 < JJ 



(4.16) 
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where 



i.e. 



fi(x) = x 



— ^,N-K+i-l e -NV(x) 



K 

n 



Hl—x 



£; — x 



x l 1 g{x) 

nf(Q-x) 



K 



g(x) 



- j-N-K -NV(x) 



I - x 



(4.17) 



(4.18) 



Now we simplify det [fi(xj)]. In order to do this we notice that from (3.8) follows that, 

K 



i tt g(x) _ ,_,k 

ei — x ' -r — 



K l-l 
el 



9{x) 



and thus we can write 

det [fi(x j )] 1<Li<K = (-) K2 det 



g(xj) 



= {-) K det det 



g(xj) 



1 



l<i,j<K 



(4.19) 



J l<j,v<K 



(-) K2 (-)^- 1) ^ i ydet 



l<i,I/<i<" 



Therefore, 



irA'/- -a _ / \K+i r m z n-k nf ef N 



x / dx\ . . . <1xk det 



g(xj) 



X j £ j, 



(4.20) 



det [W^,jv(^,^)]i< 4J <k 

i ... >< 

Now it is not difficult to calculate the last integral. Let us rewrite determinants as sums over 
permutations 



dx\ . . . dxK det 

E H*- + "' / 



g(xj) 



L 3 ci, 



dx\ 



a, p 6 Sj. 



X ... X / efa: 



l<i,j<if 

g(xi) 

xi - e CT (i) 



det [W>,jv(Mi>Zj)]i<i,.,-<K 
W/.ivC/Mi)'^) 



(4.21) 



w,' Wi,n{h p {k), x k) 

_XK-e<T{K)_ 

We co mpute the integrals above using orthogonality of monic polynomials (for details see equation 
( 4.27 ), where a similar integral is calculated). This yields 



dx\ . . . dxK det 



g(xj) 



det [Wi,N(pi,Xj)] 



l<i,j<K 



l<i,j<K 



E (-) 

a, p £ Sj. 



v„+u Pc N-K 
e a(l) 



K 



x ... x e 



N-K 
<j(K) 



n (w - Mi)) 
i 

n (w - 



(27Ti) A ^! 



is" 

n 



N-K 



K 



n ~ e ^ 



2?ri M / /7^(e<r(i),Mp(i)) 
det[W r //,jv(e 4 ,/ij)]i< 4J -< K 



(4.22) 
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We insert the obtained expression to ( 4.20 ) and with some simple algebra prove the proposition. 

□ 



IK 

4.3 Correlation function Tf u (w,GS) = [U z N [^,H\z N [ Uj ,H\ 



H 

Similar to the previous cases, the correlation functions that contain an even number of charac- 
teristic polynomials in the denominator are governed by a (different) two-point kernel. However, 
such correlation functions are not readily expressible as a determinant of a kernel divided by two 
Vandermonde determinants. Now an exact formula will be slightly more complicated. In fact, 
we prove that the correlation function J 7 j < II (m,uj) is a sum over permutations. 

Proposition 4.3. Define 2K dimensional vector e, 

e=(m,u>), Sme^O (4.23) 

Then the correlation function which contains an even number of characteristic polynomials in 
the denominator, = J-fj^w^ui), can be expressed as the following sum over permutations 



(4.24) 



= (-y 



K [lN-l\ 



i2if 



{2K)\ 



E 



det [Wui,N-K(£n(i), e Tr(K+j))] 1 



<i,j<K 



n s A(e ff (i), . . . , e^(x))A(e 7r (x +1 ), . . . , e^K)) 

The two-point kernel Wiu,n-k(^,oj) is constructed from the Cauchy transforms of monic or- 
thogonal polynomials, 

h N - K (e)h N - K -i{u) - h N _ K _ 1 (e)h N _ K (uj) 



W/7/,jv-Jf(e,w) = 
The constant jn-i is determined by equation (p? 



(4.25) 



Proof. We f ollow the procedure applied previously to the correlation function Tf\{t, fi). Instead 
of equation ( I.2C ) we obtain 



^ni(^) — i~) K Cn-2K,k 



AT! 



J N — 2K 



K 2K-N 



m4 



(N - 2K)\(2K)\ Z N A(e) 



x / dx\ . . . dx2K 

r N-2K p -NV(xi) 

x det 



K 2K 

n n ( x s - x ti 

s=l l=K+l 



X\ 64 



det \Wi,N-K{xi,x K+ j)] x<i , <K 



J l<i,j<2K 

To compute the integral above we proceed as follows 

' K 2K 

dx\...dx2K Yl IT ( x s ~ x i) 

_s=ll=K+l 



(4.26) 



x det 



x N-2K e -NV(xi) 



Xi €j 



det \Wi,N-K{xi,x K+ j)} x<i , <K 



J l<i,j<2K 



E (-r- E (-r- 



7T e s 2 



u £ St 



x / dxi . . . dxK 



r N-2K e -NV(x 1 ) 
£tt(1) - Xl 



c N-2K e -NV( XK ) 
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/ 



x / dy x 



y N — 2K 



e -NV(y l)n K {xi _ yi) 



^n(K) - J/1 



X ... X 



x / dy 



K 



^(2K) - VK 



W ItN - K {x a ( K ),y K ) 



We have 



/ 



dy 



yN-2K e -NV ( y) U K ixi _ y) 

ei-y 



Wi >N - K (xj,y) 



I 



,N-2K„-NV(y) 



dy 



lit 1 (*i -v) W + i(xi-y) 



ti-y 



X \-KN-K{Xj)'K N -K-l{y) - TV N - K (y)lV N -K-l(xj) 



3-1 



K 



3+1 



x / dy 



y N-2K e -NV(y) 



(4.27) 



\^N-K{Xj)lTN-K-l{y) - 7TN-K{y)^N-K-l{Xj) 



where the orthogonality of the monic polynomials with respect to the weight function exp [— NV(x)] 
is used. Therefore, 



y N-2K e -NV { y )U K {xi _ y) 



a-y 



K 



Wi iN - K (xj,y) 



Y[(xi - Ci) 



This yields 

J dxi . . . dx 2 x 

' x N-2K e -NV{xi) 



Wii iN - K (ei,Xj) 



K 2K 

n n i x '- x i) 

s=ll=K+l 



x det 



Xi — €j 



det [W I , N -K(x t ,x K +j)} 1 < ij < K 



J l<i,j<2K 



(-2ni) K J2 (-)"" 

rr e S 2 k 



' K 

n 



N—2K 



E (-)"• 

J a e Sk 



dx, 



(1)- 



N-2K -NV(x a(1) ) 



111 OMl) - e 7r(K+l)) 



£ 7r(o-(l)) - ^cr(l) 



W / //,Ar^K(e T (K+<j(l)),2 ; ( T(l)) 



X ... X 



/ 



dX r 



x n 2K e _ NV ^ (K)) r nf {Xa{K) _ e<K+l)) 



(K)- 



WlI i N-K(^Tr(K+a(K)),X (7 ^ K )) 
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irES 2 



K 



n N-2K 
e 7r(K+l 



E (-)"• 

- a G 5 K 



K 

n 



N-2K 



A 



I| (^(l)) - ^(K+l)) 



Wni^N-K{^Tr(K+a(l))^TT(cr(l))) 



X ... X 



A" 



n ( e *(°{K)) - ^{K+i)) 



WiiIM-K{f-TT{K+<j(K))^TT(<j(K))) 



2 A' 



(-2?ri) 



E (- 



2 A' 

n 



N-2K 



A 



17 (e ff (») - £x(A+j)) 
»>i=i 



(-2?ri) 



2 A 



2 A 

n 



N—2K 



E 



J 7T e s 2 



det [Wiu t N-K(en(i),e„(R + j))] 1 < i)J -< K 

det [W / 777 i Ar_K(e 7r ( i) , e 7r(K+J ))] ^^^^ 
A( 



We insert the above expression to the formula (4.26) and prove the proposition 



□ 



4-4 Correlation function T fy(c, p,) 



n znWh) 

~M 



, K > M 



H 



As a result of the fact that the numbers of characteristic polynomials in the denominator and 
the numerator are not equal to each other, the correlation function of characteristic polynomials 
turns out to be determined by two kernel functions. 

Proposition 4.4. Let Qrn 6j ^ 0. For the correlation function which contains K characteristic 
polynomials in the numerator and M < K characteristic polynomials in the denominator the 
following formula holds 



^(e,A) = H = 



nN-M+L-X 
N 11 



(4.28) 



A(ei, . . . , cm] hl-m+i, ■ ■ ■ > Ma) 

A 2 (e) A 2 (/7i_ M +i, • ■ ■ , hk)A(hi, Hl-m) 



x det 



Wzi,iV-M+i( e l) MX-M+l) 

Wu,n-m+l(£m, Ml-m+i) 

Wl t N—M+L (Ml i ML-M+l) 



N — M+L 



(£i,Ma) 



Wfj,jv-Af+i( c M) Ma) 

Wi-./V-M+i^l^A-) 



W7,jv-M+l(ML-M, fJ-L-M+l) ■■■ Wl ,N -M+mMX-M , Ma) 

where 2L = K + M fie. £/ie £otaZ num ber of char acteristic polynomials is even) and the kernel 
functions are defined by equations and (^AC). 
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The proof of the proposition above can be given by the same method as for the correlation 
functions !F^(e,fi) and T^ u (tz, ui). 



. 5 Correlation function Ty (e, fi) = 



n z N \pi,H] 



K < M 



H 



We also have found the representation for the correlation function Ty (e, fi) in terms of the 
kernels. Similar to the correlation function Tfy(e, jl) the correlation function Ty (e, p.) is de- 
termined by two kernels. Both these kernels now include Cauchy transforms of the orthogonal 
polynomials. Here we present the formula without a proof. 



Proposition 4.5. Let Ssm ej ^ 0. Then the following formula holds 



^(e,A) = H 



E 



b L N- 



M+L-lJ 



1 1 



nN-1 
N—M+L 7s 



A 2 (^) Ml 

A (jli, . . . , Hk\ £^(1 | M - K ) - ■■■■> e 7r(M) 



(1)' 



x det 



Wttt ,iV-M+Z, (^(l) > e ,r(l | M ~ K )) 



WlI,N-M+L(e w (M), Ml) 

Wzi,iV-M+i(e7r(M))/ 1 K') 
W777,W-M+i(e7r(l)> e 7r(M)) 



(4.29) 



Wttt M+r,(e^( M-K y £^(1 1 M ~ K )) ■•• ^JJJ,iV-M+i/( e ff ( M^y ); £ tt(M)) 

where 2L = K + M ('i.e. f/ie ioiaZ num ber of the c harac teristic polynomials is even) and the 
kernel functions are defined by equations and ( 4-25 ) 



. 6 Formula for the general correlation function 



It is possible to derive an exact expression which is valid for all the five cases considered above. 
Namely, in [pTj we have proved that 



K 

i n z n [m,H} x 

% — 

3 = 1 H 



oiV-1 
llj=N-M Ij 



x det 



A(A)A(e) 

/liV-M(ei) h,N-M+l{ci) 
7TjV-JW"(Atl) TiV-M+l(Ml) 



hN+K-l(tM) 



(4.30) 



TTJV-m(Mk) TTJV-Af+I (M-K") ■■■ ^N+K-\{^k) 

Here X and M are arbitrary positive integers. Thus formula ( 4.30| ) is also valid when the total 



number of characteristic polynomial is odd. However, the new formulae obtained in this section 
reveal a kernel structure which makes them more convenient for investigations of the asymptotic 
behaviour. 
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5. Correlation Functions of Characteristic Polynomials and the Riemann-Hilbert 

Problem 



In this section we establish a relation between the correlation functions of characteristic poly- 
nomials and the Riemann-Hilbert problem. This relation is crucial as it enables us to study the 
asymptotic behaviour for the non-Gaussian case. Below we review some aspects of the steepest- 
descent/stationary phase method for Riemann-Hilbert problems introduced by Deift and Zhou 
(a detailed presentation can be found in the book by Deift 

The method is then applied in section 6 for extracting the asymptotics of the kernels Wi t N+K(\ A*) , 
Wn,N ( £ j A*)) an d Wiii^n+k^, w )- These results combined with the propositions proved in the 
previous section give us access to the asymptotic values for all five correlation functions of char- 
acteristic polynomials (2.12|)-(pTT7). 



5.1 Relation to the Riemann-Hilbert problem. 



A technique of integrable systems (called the Riemann-Hilbert problem technique) is exploited 
in a large number of problems in mathematics and mathematical physics (see, for example, pOfl). 
The works of Fokas, Its and Kitaev ^2| relate orthogonal polynomials and their Cauchy 
transforms with an appropriate 2x2 matrix Riemann-Hilbert problem thus opening a possibility 
to apply the Riemann-Hilbert techniques to orthogonal polynomials and to the theory of random 
matrices. The observation of Fokas, Its and Kitaev (combined with semi-classical methods for 
the analysis of Riemann-Hilbert problems ||, p2[-|l8|) enabled one to understand the semi- 
classical asymptotics of the orthogonal polynomials and to provide an elegant proof of the Dyson 
universality conjecture for the Random Matrix Theory - fig)] , ||. 

The exact formulas obtained in the previous sections show that only monic orthogonal poly- 
nomials and their Cauchy transforms enter various correlation functions of characteristic polyno- 
mials. For this reason the Riemann-Hilbert problem technique arise quite naturally in the study 
of correlation functions of that type. 

All asymptotic questions we are going to address turn out to be tractable in the framework 
of the Riemann-Hilbert problem for orthogonal polynomials (Fokas, Its and Kitaev |^l|, |22j). 
Assume that the contour ^2 = M is oriented from left to right. The upper side of the complex 
plane with respect to the contour will be called the positive side and the lower part - the negative 
side. Once the integer n > is fixed, the Riemann-Hilbert problem is to find a 2 x 2 matrix 
valued function Y = Y^ (z) such that the following conditions are satisfied 



K '(z)=Yl \z) , ze 



z" 




as z 



y(")(z) - analitic in C/M 
1 

v 

• Y»(z) i > (I + Oiz- 1 )) ( 
Here Y±~\z) denotes the limit of yM (z') as z' z G R from the positive/negative side. It was 



proved by Fokas, Its and Kitaev |21j, |22j that the solution of this Riemann-Hilbert problem is 
unique and is expressed as 

Y {n \z)=( n n{z) h„(z) \ z£C/R (51) 

V ln-lK n -l{Z) Jn-lh n -l(Z) J 

where n n (z) and h n (z) are monic orthogonal polynomials and their Cauchy transforms, respec- 
tively, and the constant 7„_i is defined by equation (pT§|). 

From the propositions of the previous section we immediately observe that the correlation 
functions of characteristic polynomials are expressible in terms of the elements ofY^ n \z). More- 
over, the constant 7jv_i which determines pre-factors in the exact expressions for the correlation 
functions (see the propositions (4.1)-(4.5) and the formula ( 4.30| )) emerges in the solution of the 
Riemann-Hilbert problem as well. 
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To demonstrate the relation between the correlation functions of characteristic polynomials 
and the Riemann-Hilbert problem proposed by Fokas, Its and Kitaev it is instructive to consider 
the following example. Let us define the function 



/ detQx-g) 
^ l j ~\dct(e-tf) 



z e 



(5.2) 



H 7 e= /!= z 

The function T(z) is analytic in the whole complex plane, and J~{z) = 1, Vz G C. On one hand, 
proposition (4.2) implies that 

ir N (z) h N (z) 
In-i^n-i(z) lN-\h N -\{z) 



F{z) = dct 



(5.3) 



= det[y (JV) (z)] = 1, VzeC 

where y W (z) solves the Riemann-Hilbert problem for the orthogonal polynomials defined above. 
On the other hand, the fact that det[F (JV )(z)] = 1, Vz e C can be directly obtained from the 
definition of the Riemann-Hilbert problem (Blcher and Its §). Indeed, if Y^ N \z) solves the 
Riemann-Hilbert problem for the orthogonal polynomials it should be that 



det y w 



(*) = 



det y w 



(z) det 



l e -NV(z) 

1 



(5.4) 



det yw 



Therefore, det [Y<- N ^(z)] is analytic in C and det [Y^ N \z)] = 1 + 0(z- r ) as 
have det [Y^(z)] =1. 



Then we 



5.2 The Deift-Zhou deformations of Riemann-Hilbert problems. 



A key ingredient of the Deift-Zhou approach to extracting the asymptotics of the Riemann- 
Hilbert problem is the notion of the equilibrium measure (see II, @ H). The equilibrium 
measure is the solution of the following energy minimization problem. Assume that the value 
E is defined by 



V{s)dii{s) 



log|s — t\ 1 dfi(s)dfi(t) 



(5.5) 



Then the energy minimization problem is to find a measure d/i(s) which minimizes E . On the 
real line the equilibrium measure dfj,(x) = ip{x)dx can be uniquely determined and satisfies the 
following Euler-Lagrange variational conditions: 
There exists a real constant I such that 



log \x — y\ l dfi(y) + V(x) >l, Vie 



(5.6) 



2 / log \x - y^dfiiy) + V(x) = I, ^{x) > 



(5.7) 

For the potential V(x) — x 2m the solution of the energy minimization problem described above 
is given by the equations: 



dfi(x) = ip(x)dx, tp(x) = —{x 2 - a 2 ) 1 / 2 h 1 (x)x(- a ,a) 



hi(x) =x 2m ~ 2 + J2 X 



2m-2-2j„ 2j 



3 = 1 



n 

i=i 



21-1 



(5.8) 
(5.9) 



a = m 



n 



21 - 1 
21 



-1/2 



where X(a.~a) is the characteristic function of the interval (—a, a). 
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When the equilibrium measure is found we can define the following function : 

g(z) = J log (2: — s)ip(s)ds, s£l, Sm z =/= 

Here we take the principal branch of the logarithm, i.e. 
log(z — s) = log \z — s\+i arg(z — s) 



where 



< arg(z - s) < 7T, set, 3m z > 



(5.10) 

(5.11) 
(5.12) 



-7T < arg(z - s) < 0, set, 3m z < 

The function g(z) has the following analytical properties 
• g(z) is analytic in C \ (— oo, a] 



(5.13) 



• 9±( z ) — J log l z — s|^(s)rfs i z < —a 

a 

• 9±( z ) = J l°g l z — s|^i>(s)ds ± «7r / ip(s)ds, —a < z < a 

z 

• <?(z) = log z + C(z _1 ), z — > 00, z G C \ (— 00, a] 

Once the function <?(z) and the constant I are specified, the following transformation is intro- 



duced: 



rW(z) = e-^m (1 \z)e^e N ^ 



m$ {z)e N ^ m$ {z)e~ Nl ~ N ^ 
m$(z)e Nl + N ^ mg ) (z)e- N 3( z ) 



(5.14) 



We note that e Ng ^ is analytic in C \ [-a, a] and e Ng ^ = z N '(1 + as z -> 00. It then 

follows that the matrix function m^ 1 \z) satisfies the conditions 

• m { - 1 \z) - analytic in C/M 

fl) (1) ( e N{g-{z)-g+{z)) £ N (g+(z)+g-(z)+l-V (z)) 

• m + '(z) = m_'(z) 







,N(g+(z)-g-(z)) 



• m { - 1 \z) =1+ Oiz- 1 ) as z -> 00 

Now the analytical properties of the function g(z), the Eulcr-Lagrange variational conditions, 
and the explicit form of the function tp(s) can be exploited altogether to derive the following 
representation for the jump matrix (z) of the Riemann-Hilbert problem above: 



v^(z) 



-2mN f (t 2 -a 2 ) 1/2 hi(t)dt 

1 e la 
1 

/ -2mAT/(t 2 -a 2 ) I f /2 hi(t)(it 

e 



1 e 




e 

-2mN f (t 2 -a 2 ) 1/2 hi(t)tit 

— a 

1 



2mN 5{t 2 -a 2 Y/ 2 h 1 {t)dt 



z < —a 



z G [—a, a] 



z > a 



(In the above formula (t 2 — a 2 ) 1 / 2 has positive (negative) value when z > a (z < a).) Let us set 

z 

(p(z) =m f(t 2 - a 2 ) 1/2 /ii(t)^, z G C \ [-a, a] (5.15) 



The function tp(z) is not well-defined as it depends on the path of integration, but e ±Ntp ^ is 
well-defined and analytic in C \ [—a, a]. With help of the function <^(z) we obtain the following 
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factorization of the jump matrix (z): 

' ( ! e -2N V (z) 

1 

1 W 1 W 1 
e 2N^(z) l Jy-i J \ e 2N i>+{*) 1 

-2mN J (t 2 -a 2 ) 1/2 /ii(t)dt 



vW(z) 



1 e 




z < —a 



z E [-a, a] 



z > a 



1 



Once the jump matrix (z) is factorized we can define the new matrix- valued function 
as shown in Figure 1. 



m M(z) =m^(z) 



_ e 2N V+ (z) 1 



m^{z)=m^{z)( Ke2N l_ (z) \ 




( 2 )(z) = mF>{z) 



Fig. 1 Definition of m (2) (z) 




Fig. 2 The R-H problem for m (2) (z) 
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The matrix to*- 2 -* (z) is the solution of the Riemann-Hilbert problem on the extended contour E2 
(see Fig. 2): 

• mP^{z) is analytic in C/E2 

• (z) = rrS 2 ^ (z)v^ (z) 

• m^{z) — > / as z — > 00 

The solution (z) of the Riemann-Hilbert problem defined above has the following property: 

Proposition 5.1. Let x S (—a, a). Then for z in the vicinity of x the solution of the Riemann- 
Hilbert problem wS 2 \z) and its derivative dm^(z)/dz are bounded as N — > 00. 



Once the proposition (5.1) is proved one can observe that Kjy(x, x) defined by the equation (2.7) 
is equal to Nip(x) to the leading order, i.e. 



K N (x,x) = Nip(x) + 0(1) as N 



(5.16) 



6. Asymptotics of the Kernels 



In this section we use the results outlined above to determine the asymptotic behaviour of three 
kernels Wi t N+K{\ p), Wh^n{\ A*), an d Wiii,n-k(\ m) m the Dyson's limit. This is achieved 
by three subsequent transformations. The first step is to express the kernels Wi n+k{^ AO — 
Whi,n—k(\ m) m terms of matrix elements of (z). We then rewrite them in terms of mP-'fe) 
and, finally, in terms of rrS 2 \z) defined by Fig. 1 and Fig. 2. The reason for these transformations 
is that m^ 2 '(z) and its derivative are bounded matrix valued functions as N — > 00. It is this fact 
and the equation det (z) = 1 that enable us to find the large N asymptotic of kernels in the 
Dyson's scaling limit. The obtained asymptotic formulae are summarized in Table 3. 



Kernel 


Large N asymptotic 


Wi,n+k(x,(,v) 


c% [Np(x)\ e^^e^XC+^S^C - v ) 


Wji, N (x,C,ri), 


-3L[Np(x)] e- a ^-^S„{C-v) 







Table 3. Asymptotic of kernels 



6.1 Large N limit of Wi^n+k{\ A*) 



We fix a point x G (—a, a). The interval {—a, a) is the support of the equilibrium measure for 
the given potential function V(x). For simplicity we assume that the support of the equilibrium 
measure includes only one interval. This is the case for the potential function V(x) — x 2m , 
m > 1. Introduce new coordinates £,77 

X = X+ Np\xj> » = X+ NpJx) (0) 
In what follows we consider the Dyson's scaling limit. In such a limit the difference between points 
A and p goes to zero, the size N goes to infinity, the product N(X — p) remains finite. In equation 
([O]) p(x) = Kn(x 1 x)/N is the density of states, with the kernel Kn(x,x) being given by the 
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equation ( |2.7|) . In terms of the new coordinates the kernel Wi,n+k{\ A*) = Wi,n+k(x, C v) is 
expressed as a determinant of monic orthogonal polynomials, 

Wi,if + K(x,c,n) = [N P (x)] {(-vr 1 



x det 



(6.2) 



Thus the problem about asymptotics of the two-point kernel function Wt,n+k(\ m) is reduced 
to the investigation of the large N asymptotics of the determinant in the equation (|6.2|) . 



Proposition 6.1. (Large N asymptotics of the kernel Wi t N+K{x,C->v))- 

Let x S (—a, a), —9 < £, r\ < 9, £ ^ ry and AT 3> if. TTien in the large N limit the kernel function 
W/,iV+K (x, C, Tj) is related with the universal kernel Si{(, — rf) as 

W ItN+K (xX,v) = [cn] 2K [Np(x)} e NV ^e a ^+^ [5/(C - v) + 0(1/ N)] 
where the universal kernel Si (( — rf) is 
sin(7r(C - rf)) 



tt(C - rf) 



(6.3) 
(6.4) 



Proof. Equation (5.1) enables us to write the kernel W/ jv+j<-(x, C, rf) in the form: 



Wj 



N+K 



(x, C, Tj) = [7. 



N+K-H 



[Np(x)] [C-r,Y 



x dot 



Y [N+K) (x + ^— 

+ 



'/ 



JVp(x) - 



(6.5) 



The large TV asymptotics is completely determined by the determinant^ in the equation above. 
Large N limit of this determinant is considered in Deift's book [|l2| , Chapter 8. Here we reproduce 
his derivation. Equation flS.14 ) gives 



det 



Y^ +K \x 



Y 



Np(x 



Y} N+K \x 



Np(x 



Np(x 

N[g + (x+(/Np(x))+g + (x+r,/Np(x))+l] 



Np(x) 



(6.6) 



x det 



m 21 ( X + Np(x) ) m 21 ( X + Np\x) ) 

where mS^[z) is the solution of the transformed Riemann-Hilbcrt problem (see the previous 
section). Rewrite the determinant in terms of rrv- 2 >(z) (see Fig. 1 and Fig. 2 where m^(z) is 
related with mP^(z)). This leads to the following expression for the kernel 

wi, N+K (x, c, rf) = [jN+K-ir 1 [Np{x)} [c - vr 1 



x e 



N[g+{x+C/N p{x))+g+(x+r,/N p(x))+l] 



(6.7) 







e 2N v+ (xr,) 






+ 










+ 



x det 

where x c = x + .t„ = x + 

matrix valued functions in the vicinity of the point x and det m^ 2 \z) = 1. This observation 



Remember that m^(z) and its derivative are bounded 



^Since N 2> K, we do not distinguish between N and N + K when studying the asymptotics. 
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enables one to expand all functions in the determinant in the vicinity of the point x, and to find 
the leading term: 

w ItN+K ( X , c, v) = [in+k-iV 1 [n P {x)\ [c - vr 1 



x e 



N [g+ (x+C/Np(x))+g + (x+ V / N p(x))+l] 



(6.8) 



e 2N V+ (x+ri/Np(x)) _ e 2N V+ (x+(/Np(x)) 



[i + o^- 1 )] 



The next step is to express the functions g+(x + Q/N p(x)), g+{x + rj/Np(x)) in terms of the 
functions tp+(x£), ip+(x. q ). It can be done as follows. Let z G (—a, a). When z G (—a, a) 
the following equations hold (see properties of the function g±{z) summarized in the previous 
section): 



9+{z) + g-{z) 



log \z — s\Tp(s)ds 



(6.9) 



9+{z) ~ 9-i z ) = 2«7r / tfj(s)ds 



(6.10) 



We use the Euler-Lagrange variational condition ( |5.7| ) together with the equations above and 
find that the functions g±(z) are completely determined by the function tp+(z), by the potential 
V(z) and the constant I. Namely, 



(g+(z) = -<p + (z) + V(z)-± 
\g-(z)= V+ {z) + V{z)-± 

This gives 

W l!N+K (x,C,r,) = hN+K-iV 1 [Np(x)] [C - vV 1 e NV ^e a ^+^ 



(6.11) 



,N[ V+ (x+C/Np(xy)- V+ (x+r,/Np(x))] 



In the vicinity of the point x we have 



-N[ v+ (x+ v /Np(x))~ v+ (x+(/Np{x))] 



(p+(x + 



c 



-) = (f+(x) + iirtp(x)- 



Np(x) 'Np(x 

A similar expression is obtained for f+(x J - 

( 5.16| ) and the relation Kn(x,x) = Np(x) 



Np(x 



O 



N 2 



(6.12) 



(6.13) 



We complete the proof using the equation 

□ 



6.2 Large N limit ofWu^i^,^) 



In this subsection we investigate Dyson's scaling limit of the kernel W/_r,Ar(e, p) constructed from 
monic orthogonal polynomials and their Cauchy transforms. Similar to the procedure applied to 
the kernel function Wi^n+k{^, p) we introduce new coordinates, 

, c 



p. = x + 



V 



Np{x) ' n Np{x) 
In the coordinates £, r\ the kernel Wh,n (e, p) = Wh : n{x 1 (, rj) has the form: 



(6.14) 



x det 



(6.15) 



We can see that in order to derive an asymptotic expression for the kernel Wii_n(x, £, v) we 
need the asymptotics of the determinant which contains monic orthogonal polynomials and their 
Cauchy transform. The Riemann-Hilbert technique proves to be a convenient tool in this case 
as well. 
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Proposition 6.2. (Large N asymptotics for the kernel Wn t N(x,CiV)) 

Let x G (—a, a), —9 < £, rj < 6, C V> ^ m C 0. Then the following asymptotic expression for 
the kernel Wji t N(x,^,rj) holds: 



W n , N {x, £ V) = [Np(x)} e - Q (^-") [S n (<: - rf) + 0(1/ N)] 

The universal two-point kernel S7/(C — rf) is expressed as 



c-n 



C-v 



Sto £ > 0, 
3m(<0 



(6.16) 



(6.17) 



Proof. Express the determinant in the equation (6.15) in terms of Y^ N \z) which is the matrix 
valued solution of the Riemann Hilbert problem for the orthogonal polynomials. We have 

W n , N (x, C, rj) = [Np(x)] [r, - CP 1 fr^-ip 1 



x det 



Y ll ( X + Np(x) ) Y 12 ( X + Np(x) ) 

i C ) 



Y 2 f\x- 



) Y$\x + 



Np(x) ■ 



, , (6-18) 

Np(x)> ± 22 \* ~r Np{x)' 

Equation ( 5.14 ) enables one to replace the elements of Y^ N ' by the elements of mS 1 '(z) (we 
remind that mP''(z) is the solution of the transformed Riemann-Hilbcrt problem). Then we 
obtain 

Wn, N {x, C, rf) = [Np{x)} [n - fl -1 frjv-i] -1 e N ^ x +^ N P( x ))-9( x +C/Np(x))] 



x det 



(6.19) 



TO 21 ( X + Np(x) ) TO 22 ( x + Np(x) ) 

The elements of m^(z) can be further replaced by the elements of mS 2 \z) (see Fig. 1 and Fig. 
2 where the Riemann-Hilbert problem for (z) is specified and the relation between mW (z) 
and mS 2 \z) is shown). Again, m^ 2 \z) and its derivative are bounded near the point x, and 
det m^(z) = 1. Thus we can expand around the point x and show that the determinant in 
equation ( 6.19| ) is equal to one to the leading order. Therefore, 

W n , N {x, C, rj) = [Np{x)\ [ v - CP* llN-i]' 1 e N M*+v/Kp(*))-9(*+t/"p(*))] [i + 0(1/N)} 

Introduce z\ = x + r)/Np(x), Z2 = x + Q/N p(x). Then we have 



N[g(z 1 )-g{z 2 )] =N 



log(zi — s)ip(s)ds — / log(z2 — s)ip(s)ds 



N / lo. 



Zl - s 



N J log 1 

— rr lim 

p(x) A^O 

[ri-Q 



z 2 - s 



Z 2 - S 



ip(s)ds 

Tp(s)ds 



ijj(s)ds 
: — s ±i5 



p(x) 



irHip(x) =F mij}(xj\ 



(6.20) 



0(1/N) 



where + (— ) corresponds to the positive (negative) imaginary part of £, and Hip(x) stands for 
the Hilbert transform of ip(x), 



Hi/j(x) = - P.V. 

7T 



ip(s)ds 



(6.21) 
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Furthermore, from equation (|5.7| ) (the second Euler-Lagrangc condition) we observe that 

HiP(x) = ± V'(x) (6.22) 
Now use that ip(x) is equal to the density of states in the large N limit to obtain 

N [g( Zl ) - g(z 2 )} = ±i7r(C - r,) ^||(C ~v) + 0(±) (6-23) 

which completes the proof. □ 
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6.3 Large N limit ofWiu^N-K^,^) 



Finally we investigate the Dyson's scaling limit for the kernel Wjti,n-k(^> u). This kernel is 
constructed from the Cauchy transforms of monic orthogonal polynomials, and the transforms are 
not analytic. For this reason we need to consider different situations corresponding to different 
signs of the imaginary parts. Just as before we introduce new coordinates 

e = x + C/Np(x), LU = x + rj/Np(x) (6.24) 

where 3m ( ^ 0, 3m rj ^ 0. In terms of these coordinates we define Wjjj t jf— k(x, Ci if) = 
W IU m-k{x + j^-y x + 



Proposition 6.3. (Large N asymptotics of the kernel Wih.n-k (x.X,v)) ■ Letx E (—a, a), —0 < 
C, ry < 6, 3m C, ^ 0, 3m rj ^ 0. Then the large N limit of the kernel function Wui,n-k(x, C v) 
is related to a universal kernel Sni(£,rf) as 

Wiii,n-k(x,C,ii) 



= M" 1 [Np(x)\ e -wv(») e -a(*)(c+,) [ 5jJj(Cj v) + oiN- 1 )] (6.25) 
where the universal kernel Sm (£, rj) is given by the formula 

3m £ > 0, 3m n < 0, 

S, i ,(,.//) { ^ 3m C < 0, 3m V > 0, (6.26) 
otherwise 



Proof. We give the proof only for the case 3m £ > 0, 3m r/ < 0. The other three cases 
with different signs of imaginary parts can be considered in a similar manner. The kernel 
Wiii,n-k{x,CiV) can be expressed as 

Win, N - K (x, C, v) = [N P (x)] [C - vV 1 



xdet h N - K {x + (/Np{x)) h N „ K {x + n/Np{x)) 
Hn-k-i(x + C/ n p( x )) h N - K -i(x + i]/Np(x)) 

Now we exploit the relation to the Riemann-Hilbert problem for the orthogonal polynomials. We 
replace the Cauchy transforms in the determinant above by the corresponding elements of the 
matrix Y^ N ~ K > which is the solution of the Riemann-Hilbert problem. We have 

W iu ,n-k{x, C, V ) = [Np(x)} [C - vV 1 [in-kV 1 



x det 



Y^- K \x + d/Np{x)) 
Y^- K \x + C/Np(x)) 



Y} 2 N - K \x + V /Np(x)) 
Y^- K \x + V /Np(x)) 



Then we employ the transformation (equation ( 5.14 )) from the solution Y^ N ~ K '(z) of the original 
Riemann-Hilbert problem to that of the new Riemann-Hilbert problem mS 1 '{z) defined by the 
jump matrix v^'(z) (N should be replaced by N — K). This yields 

WlII,N-K (X, C; V) = l N P(x)] [C - [IN-kY 1 



x exp [-{N -K)[l + g+(x + (/Np(x)) + g+{x + n/Np{x))]] 



x det 



rr$(x + t/Np{x)) 
m£(x + (/Np(x)) 



rrii2 (x + r)/Np(x)) 
m$(x + r)/Np(x)) 
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In turn, the function m^'(z) is related to the matrix valued function m^{z) (see Fig. 2) which 
is the solution of the deformed Ricmann-Hilbert problem defined by Fig. 2. Correspondingly, 
we rewrite the kernel in terms of the elements of m^(z), 

W m>N - K (x,C,v) = [Np(x)] [C-vV 1 [in-kY 1 



x exp [—(N -K)[l + g+{x + (/Np(x)) + g+(x + r]/Np{x))]} 



x dot 



rr$(x + S/Np(x)) 
m$(x + C/N P (x)) 



(2) 

m-L2 (x + 7]/Np(x)) 

(1) 

my(x + r]/Np(x)) 



As m^ 2 \z) and dm^ 2 ' (z) / dz are bounded matrix valued fun ction s in the vicinity of the point x, 
and g+(z) + g+(z) + I = V(z) for z 6 (—a, a) (see equation ( 6.11 )) we obtain the following large 
N limit for the kernel Win t N—K(%>GiW) 

W m ,N-K(x,C,v) = [Np(x)] [C-vV 1 [In-kY 1 e - NV ^e~ a ^<+v) 









( 2 )c \~ 




( 2 )/ ^" 






X 


det 




my (x) 


+ 


my(x) 




+ (D(l/N) 








+ 







However, in the vicinity of the point x we have (see Fig. 2) 





( 2 ) 1 \ 


+ 


(2) 
m 12 


(x)_ 




(2) 
m ll 


(x) 




(2) 
m 12 


(x) 


( 


m 21 '{x) 


+ 


(2) 

my 


(x)_ 


:)-( 


(2) 
m 21 


(x)_ 




(2) 

my 


(x) 



i.e. 



m 12W 



"Hi O) 



ra 21» 



my [x) 



Taking into account that 

"Hi W| 
w 21 '(x) 



det 



my{x) 
( 2 )c > 

™22 W 



1 

-1 



we finally obtain 

W hi ,n-k(x, C, v) = [Np(x)} [C ~ 17]- 1 [Tiv-*]" 1 e - NV ^e- a ^+^ [1 + 0(1/JV)] 



□ 



7. Negative Moments 



In this section we derive the asymptotic expression (2.24) for the negative moments of character- 
istic polynomials M% N (S) defined by the equation (2.23). The negative moments can be obtain 
as limiting values of the correlation function Tf u (x + C,/N p(x)). The large N asymptotic for 
that function is given by the equation ( 2.15| ). We define three K dimensional vectors, C + , C~, 
and 8. The components of the vectors £ + (C — ) are pure imaginary and have positive (negative) 
imaginary parts. The vector 6 has all components equal to each other and equal to the real pa- 
rameter 5 in the definition of the negative moments (equation (2.23)). We exploit our asymptotic 
formula ( |2.15 ) for the correlation function Tf II (x + Q/N p{x)) and write the following expression 
for the negative moments: 

-KNV(x) 

Jk(S) 



(-) 



K 



[ lN ] K [Np(x)f , 



(2K)\ 



(7.1) 
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where 



Jk(S) 



lim 

C±^±i<5/2 



E 



det [§/// (Cx(i), Ctt(j+k))] i 



<i,j<K 



Ltt e s k A ^ 7r ( 1 )' ■ ■ ■ 'C7r(K))A(C 7r (x+i), ■ ■ • ,Ctt(2K)) 

Following the method by Brezin and Hikami's (|, ^ we represent the determinant of the kernel 
divided by two Vandermonde determinants as a contour integral, i.e. 



E 



det [Sjjj (Cr(i), CttO'+K))] ! 



<i,j<K 



7T e s 2 



A(Ctt(1), ■ • ■ i C7r(_R-))A(C7r(i4:+l)i • ■ ■ , Ctt(2A')) 



(7.2) 



1 



E 

£ S2RT 



n 



duidvi A(w)A(w)nf $in(ui,Vi) 



A (27Ti) 2 



where the contours of integration are chosen in such a way that all components of the vector £ 
give rise to contributions to the integral as simple poles. We have to compute the integral above 
when 

C = (it/2, . . . , it/2, -i5/2, -i5/2) = (C + , C ) (7.3) 

All those permutations tt € S2K that lead to the same vector 7r£ produce the same contribution 
to the integral. This observation permits us to rewrite the right-hand side of the equation (7.2) 

as 

K[ E/ /Il (27ri) 2 K " 

1 I! [ U i-Q<T{3)\ l V i - Ca(K+j)\ 

where a 6 S2k/^>k x Sa' should be understood as permutations exchanging elements between 

the two sets, (iS/2, . . . , id/2) and (—iS/2, . . . , —iS/2) (the dimension of each set is K). Consider 

all such permutations that replace K\ < K elements of the first set, (i8/2, . . . , iS/2), by K\ 

elements of the second set, (—iS/2, . . . , —iS/2). The number of such permutations is equal to 
2 

. A precise location of new elements in the first and the second sets does not affect 
2 



K\ 



the integral, so all 



Kl 



K 

Jk(6)=K\ 



give the same contribution. Therefore, 



Ki=0 



K\ 



{K — K{)\K\\ 



T K,Ki 



(7.4) 



where 



r K,Ki 



K 

n 



dujdvj 



A(«)A(«)nf S UI ( 



■1 (27rf) 



2 K 



Y[[(u l + 4)(v l -il)] Kl [(„,_»§)(„,+»«)] 



K-Ki 



Now we rewrite the contour integral lf^' Kl as a determinant of a kernel divided by two Van- 
dermonde determinants. It follows from explicit expressions for limiting kernels summarized in 
table 2 that E>ui(a,f3) — to the leading order, if the variables a and j3 have imaginary parts 
of different signs. Then we obtain 

jf' Kl = lim 

Si->8/2 
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We insert the explicit expression for the kernel — if) (see Table 2) to the above formula 

and find that 



T K,K\ 



(-i) K Kl 
det 



x lim 

S t ->6 



A 2 



( 5i Ski_ 

I 2 5***5 2 



det (t-xt) 

i J Ki+l<i,j<K 



A 2 



(7.5) 



2 ' ■ ' ■ ' 2 

Exploiting once again Brezin and Hikami's representation in terms of a contour integral we find 
det 

= det {a mn ) 0<m _ n< M-i ( 7 - 6 ) 



lim 



S * +S = J l<i,j<M 



where 



A2 (§x §M-\ 



f &" d n 



u + v + 6 



u—v—0 



(-) m + n f(m + n)l 



Now it is straightforward to compute det (a mn ). We obtain 



lim 

6i->6 



det 



A 



di+d i J l<i,j<M 
2 I Si S M \ 



S M2 



This equation together with the equation (7__5) yields 



tK.Ki 



L 8 

r K,Ki 



{-i) K K\ 5 



(i.e. I s ' 1 does not depend on K\). Taking into account that 



K 

E 

Ki=0 



K\ 



(K-Ki)\K]\ 



(2-g)I 

(K\y 



(7.7) 



(7.8) 



and equations (7.1), ( |7.4[ ) we finally obtain our result (2.24) for the negative moments of the 
characteristic polynomials. 



8. Two-point Correlation Function 



The present section serves for illustrating the utility of the correlation functions constructed 
with help of the characteristic polynomials. We are going to demonstrate how one can derive 
the resolvent two-point correlation function from our asymptotic result for the average values of 
ratios of the characteristic polynomials. We consider the case when a(x) = = 0. Then 

the correlation function Tjj{x + (/Np(x), x + rj/Np(x)) is universal in the Dyson scaling limit. 
As the result, the answer for the resolvent correlation function will be universal as well, as a 
by-product we provide a new proof of the universality of the two-point correlation function of 
eigenvalue densities. 

The resolvent two-point correlation function is defined as 



S 2 (E,E r ) = N~ 2 ( Tr 



1 



E-H 



Tr 



1 



E'-H 



(8.1) 



H 



where the (complex) energy E has a positive imaginary part and the energy E' has a negative 
imaginary part. To consider the scaling limit we introduce new coordinates 



E = x + rj-y/N p{x), ^sm r\\ > 



(8.2) 
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E' = x + 7 l2 /Np(x), 3mry 2 <0 (8.3) 

To connect the two-point correlation function S-z(x+rii j 'N p(x) , x+r)2/Np(x)) with the correlation 
function Tff 2 (x + £/Np(x), x + fj/Np(x)) investigated above we exploit the following identity 



Tr 



= [Np(x)} 



x + r)/Np(x) - H 
where x n = x + rj/Np(x) . We have 

S 2 {x + rn/Np(x),x + rj 2 /Np(x)) 



d 71 Z N [x v ,H] 
Z N [x n ,H] 



(8.4) 



det (x Vl — H) det (x V2 — H) 
det (x£ ± — H) det (x^ 2 — H) 



= W)Y d; um Fn\x + C/Np(x), x + f,/Np(x)) 



Equation ( [2.14| ) yields 



Cl='?l,C2=f2 



FiiCx + C/Np(x), £ + fj/Np(x)) 



s 'i7r(Ci-C2) 
Cl-C2 



J*(m-m) fa ~ d)fa - C2) _ c - i7r(m - m ) (Vi - (2)^2 - Ci) 

Vl ~ V2 Vl ~ V2 

The next step is to compute derivatives. In particular, we find 

T(V1-V2) fa ~ Cl)fa ~ C2) " 

Vl~V2 Ui=Vi,Ca=V2 
e i7r(i7i-r; 2 ) 



d 2 



m - m 



and 



-nr(m-V2) fa ~ Cg)fa ~ Cl) 



Vi - V2 



ir 2 {r 1 i-V2)e~ Mm ~ r,2) + 



-i7r(?7i-jj2) 



(8.5) 



(8.6) 



vi - m 

which gives the well-known result (see |43L p3fl ) for the resolvent two-point correlation function: 



S 2 (x + m/Np(x), x + r) 2 /Np(x)) 



[irp(x)Y 



1 - 2i 



suitt^ - 771) e'^^-^ 
[tt(v2 - Vi)? 



(8.7) 



9. Summary and Discussions 



In this paper we prove three basic statements: 1) correlation functions of characteristic poly- 
nomials are governed by two-point kernels, 2) the kernels are "integrable" in the sense of the 
definition by Its, Izergin, Korepin and Slavnov [^8)-(3^], 3) the kernels are constructed from monic 
orthogonal polynomials and their Cauchy transforms. As a consequence, it becomes quite nat- 
ural to exploit a relation to the Riemann-Hilbert problem for orthogonal polynomials proposed 
by Fokas, Its and Kitaev ||, ||. 

It is known that the simplest correlation functions, i.e. the moments of the characteristic 
polynomials, can be described in terms of non-linear differential equations (see works of Forrester 
and White |2^,|24|], and also the paper by Kanzieper Q, and by Splittorff and Verbaarschot Q). 
As for more complicated correlation functions a description in terms of differential equations is 
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unknown. While in the present paper we focus on asymptotic questions, the " integrability" of 
kernels suggests that such description should be possible. 

In this paper the discussion is restricted to the ensemble of unitary invariant Hcrmitian ma- 
trices. However, the case of compact group ensembles (circular ensembles) can be approached 
in the same way. The circular ensembles are even simpler as the representation theory of com- 
pact groups (inapplicable for Hermitian random matrices) can be exploited there. Indeed, the 
method of "dual pairs" (Nonenmacher and Zirnbauer [55|, Q) enables one to find exact for- 
mulas for correlation functions of characteristic polynomials for the circular ensembles. However, 
this method is based on an interpretation of characteristic polynomials as characters of " spinor" 
group representations. As a result, its applicability is restricted to group ensembles. The method 
proposed in this paper is based neither on the representation theory, nor on specific features of 
Hermitian matrices. The only fact which is important is that the ensemble under considerations 
is of (3 = 2 symmetry class. For this reason it is more general and can be applied both to 
Hermitian random matrices and to the group ensembles. 

Correlation functions of characteristic polynomials of random matrices for the unitary circu- 
lar ensemble and the ensemble of Hermitian matrices are sometimes represented as Toeplitz / 
Hankel determinants, correspondingly. Results for Toeplitz determinants with rational generat- 
ing functions can be found in Basor and Forrester Q. However, Hankel determinants are less 
well studied. The correlation functions investigated in the present paper can be understood as 
Hankel determinants with rational generating functions. Thus, asymptotic and exact statements 
for Hankel determinants equivalent to our results should be possible to make. 

We hope that the method proposed in this paper can be modified as to provide an access 
to other symmetry classes of invariant ensembles j3 = 1,4 (for the gaussian case some attempts 
were undertaken recently in ||). Another important goal is to apply them to ensembles of 
non-Hermitian random matrices |p6[ P8| , which are certain deformations of the invariant class. 

The most challenging problem is to investigate the conditions of universality of the discussed 
correlation functions for non- invariant non-gaussian ensembles with independent, identically dis- 
tributed (i.i.d.) entries. So far the progress was rather limited and restricted to the specific 
choice of the probability measure (see Johansson |40| ). At the same time, non-rigorous heuristic 
methods hint to a kind of universality covering also the so-called ensembles of sparse random 
matrices, see Mirlin, Fyodorov Q and Fyodorov, Sommers p7[ . 

We leave a detailed investigation of these issues for future research. 
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Appendix A. 



Let us consider the sum 
<y- 7r £ S n + m 



X 9-K(l)(lpl) ■ ■ ■97r(m)(i } m) (A.l) 

x ^cr(m+l)7r(m+l) ■ ■ ■ ^a(m+n)Tr(m+n)] 
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In order to reduce this expression to a determinant form (to a determinant of size to x m) we 

, k m ] satisfying the following condition 



proceed as follows. We fix the set of numbers [k\, 

n + m > fei > hi > . . . k m > 1 

Let us denote by <t(1), . . . , cx(to) ( and by 7r(l), 
numbers 1,2, ... , to e nd u p inside the set [fei, hi, 
the sum S (equation (A.l)) can be rewritten as 



(A.2) 

,7r(m)) the permutations under which the 
,k m ]. Once such notations are introduced, 



S = 



E 

fei>...>fc m >i 



E 

O", 7T G §n + r> 



(-) u * + *[fa(i)(cpi)---U( m ){¥>m) 



■ <?7r(m) (i>rn) 



(A.3) 



X ^5-(m+l)7r(m+l) ■ ■ ■ ^d(m+n)Tr(m+n)] 

Under permutations cr, 7r the numbers to + 1, . . . ,m + n remain outside of the set [fci, . . . , fe m ]. 
The sets [1, . . . , m], [to + 1, . . . , m + n] do not mix, and this leads to the following expression for 
S: 

n+m 



s = 



E 



fei>...>fc m >l Lct g s m 

E (-)' y *3fe*(i)(V'i)---gfc s(m) (V'm) 



(A.4) 



.7fe s„ 



7^ ( _ ) 1/cr+1/7r ^a-(m+l)7r(m+l) 
cr, 7T £ S n 



• &a(m-{-n)7v(m-\-n) 



From the above representation we immediately conclude that 

n+m 

S = n\ det ^ det ^ (^')] (A.5) 

fci>...>fe m >i 

where the indices i, j take values from 1 to m. Now the sum S can be rewritten as a determinant 
to x to, as a consequence of the formula 



^ det [/fc<(^-)] det [fffc 4 (V'j)] = dct 
fei>...>fc m >i 



n+m— 1 



E h(<Pi)9x(i>j) 



A=0 



(A.6) 



(This fact is a generalization of the theorem on determinants. This theorem ( which can be found 
in the book of Hua |B4j) states that 



det 



+ 00 



,A=0 



Edet \A\, i], „. . , det 



l<i,j<m 



Ai>...>A m >0 



l<i,j<m 



(A.7) 



A nice proof of the formula ( A.6 ) can be found, for example, in the paper by Balantekin and 
Cassak 101). 



References 

[1] Andreev, A., V., and Simons, B., D.: Correlators of spectral determinants in quantum chaos. Phys. Rev. 
Lett. 75, 2304-2307 1995 

[2] Balantekin, A., B., and Cassak, P: Character expansions for the orthogonal and symplectic groups. J. Math. 
Phys. 43 (1), 604-620 2002 

[3] Basor, E., L., and Forrester, P., J.: Formulas for the evaluation of Toeplitz determinants with rational 

generating functions. Mathematishe Nachrichtcn, 170, 5-18 1994 
[4] Berry, M., V., and Keating, J., P.: Clusters of near- degenerate levels dominate negative moments of spectral 

determinants. J. Phys. A: Math. Gen. 35, L1-L6 2002 
[5] Bleher, P., and Its, A., R.: Semiclassical asymptotics of orthogonal polynomials, Riemann- Hilbert problems, 

and universality in the matrix model. Ann. Mathematics 150, 185-266 1999 
[6] Brezin, E., and Hikami, S.: Characteristic polynomials of random matrices. Commun. Math. Phys. 214, 

111-135 2000 



34 R-H APPROACH FOR CORREL. OF CHARACTERISTIC POLYNOMIALS 

[7] Brezin, E., and Hikami, S.: Characteristic polynomials of random matrices at edge singularities, Phys. Rev. 
E 62 (3), 3558-3567 2000 

[8] Brezin, E, and Hikami. S.: New correlation functions for random matrices and integrals over super groups, 



arXiv:math-ph/020800l| 2002 
[9] Conrey, J. B., Farmer, D. W., Keating, J. P., Rubinshtein, M. O., and Snaith, N.C.: Autocorrelation of 

Random Matrix Polynomials, arXiv: mat.nt/0208007 2002 
[10] Conrey, J. B., Farmer, D. W., Keating, J. P., Rubinshtein, M. O., and Snaith, N.C.: Integral moments of 

zeta- and L-functions, arXiv: mat.nt/0206018 2002 
[11] Deift, P.: Integrable operators, Differential operators and spectral theory: M. Sh. Birman's 70th anniversary 

collection (V. Buslaev, M. Solomyak, D. Yafaev, eds.), American Mathematical Society Translation, ser.2, v. 

189, Providence 69 1999 

[12] Deift, P.: Orthogonal Polynomials and Random Matrices: A Riemann-Hilbert Approach. Courant lecture 
notes, 3 2000 

[13] Deift, P., Its, A. R. and Zhou, X. A.: A Riemann-Hilbert approach to asymptotic problems arising in the 
theory of random matrix models and also in the theory of integrable statistical mechanics. Ann. of Math. (2) 
146, 149-235 1997 

[14] Deift, P. and Zhou, X. A.: A steepest descent method for oscillatory Riemann-Hilbert problem, Asymptotics 

for the MKdV equation. Ann. of Math. (2) 137, 295-368 1993 
[15] Deift, P. and Zhou, X. A.: Asymptotics for the Painleve II equation. Comm. Pure Appl. Math. 48, 277-337 

1995 

[16] Deift, P., Venakides, S. and Zhou, X. A.: The collisionless shock region for the long-time behavior of solutions 

of the KdV equation. Comm. Pure Appl. Math. 47 (2), 199-206 1994 
[17] Deift, P., Kriecherbauer, T., McLaughlin, K. T-R., Venakides, S. and Zhou, X. A.: Strong asymptotics of 

orthogonal polynomials with respect to exponential weights. Comm. Pure Appl. Math. 52 (12), 1491-1552 

1999 

[18] Deift, P., Kriecherbauer, T., McLaughlin, K. T-R., Venakides, S. and Zhou, X. A.: Uniform asymptotics for 
polynomials orthogonal with respect to varying exponential weights and applications to universality questions 
in random matrix theory. Comm. Pure Appl. Math. 52(11), 1335-1425 1999 
[19] Efetov, K., B.: Supersymmtry in Disorder and Chaos. Cambridge University Press, Cambridge 1997 
[20] Fokas, A., S., and Zaharov, V., E.: (eds.) Important developments in soliton theory. Springer- Verlag, Berlin 
1993 

[21] Fokas, A., S., Its A., R., and Kitaev A., V.: An isomondromy approach to the theory of two-dimensional 

quantum gravity. Russian Math. Surveys 45 (6), 155-157 1990 
[22] Fokas, A., S., Its A.,R., and Kitaev A.,V.: Discrete Painleve equations and their appearance in quantum 

gravity. Comm. Math. Phys. 142 (2), 313-344 1991 
[23] Forrester, P., J., and Witte, N., S.: Application of the tau-function theory of Painleve equations to random 

matrices: PIV, PII and the GUE Comm. Math. Phys. 219 (2) 357-398, 2001 
[24] Forrester, P., J., and White, N., S.: Commun. Pure Appl. Math. 55, 679-727 2002 

[25] Fyodorov, Y., V.: Negative moments of characteristic polynomials of random matrices: Ingham-Siegel in- 
tegral as an alternative to Hubbard- Stratonovich transformation. Nucl. Phys. B [PM] 621, 643-674 2002 



arXiv:math-ph/0106006 



[26] Fyodorov, Y., V., and Khoruzhenko, B., A.: Systematic Analytic Approach to Correlation Fuctions of 

Resonances in Quantum Chaotic Scattering. Phys. Re v. Lett. 83, 65-68 1999 
[27] Fyodorov Y.V., and Sommers, H.-J. Universality of "level curvature" distribution for large random matrices: 

systematic analytical approaches Eur. Phys. J B 99 (1): 123-135 1995 
[28] Fyodorov, Y., V., a nd Sommers, H-J: Random Matrices close to Hermitian or unitary: overview of methods 
and results. arXiv: rilin.CD/020705l| 2002 



[29] Fyodorov, Y., V., and Strahov, E.: Characteristic polynomials of random Hermitian matrices and 
Duistermaat- Heckman localisation on non-compact Kahler manifolds Nucl. Phys. B [PM] 630, 453-491 2002 

[30] Fyodorov, Y., V., and Str ahov, E.: On corr elation functions of characteristic polynomials for chiral Gaussian 
Unitary Ensemble arXiv: hcp-th/0205215 2002 

[31] Fyodorov, Y., V., and Stra hov. F.: An exact formula for general spectral correlation function of random 
Hermitian matrices arXiv: math-ph / 020405 1| 2002 



[32] Gonek, S., M.: On negative moments of the Riemann zeta-fuction. Mathematika 36, 71-88 1989 
[33] Haakc, F.: Quantum Signatures of Chaos. 2nd ed., Springer 2000 

[34] Hua, L., K.: Harmonic Analysis of Functions of Several Complex Variables in Classical Domains. American 

Mathematical Society, Providence, Rhode Island 1963 
[35] Hughes, C, P., Keating, J., P., and O'Connell, N.: Random matrix theory and the derivative of the Riemann 

zeta function. Proc. R. Soc. Lond. A 456, 2611-2627 2000 
[36] Hughes, C, P., Keating, J., P., and O'Connell, N.: On the characteristic polynomials of a random unitary 

matrix. Commun. Math. Phys. 220(2), 429-451 2001 
[37] Its, A., R., Izergin, A., C, Korepin, V., E.: Themperature correlators of the impenetrable Bose gas as an 

integrable system. Comm. Math. Phys. 129, 205-222 1990 
[38] Its, A., R., Izergin, A., G., Korepin, V., E., and Slavnov, N. A. : Differential equations for quantum 

correlation functions. Int. J. Mod. Phys. B4, 1003-1037 1990 
[39] Its, A., R,., Izergin, A., C, Korepin, V., E., and Slavnov, N. A. : The quantum correlation function as the 

t function of classical differential equations. Important developments in soliton theory, A. S. Fokas and V. 

E. Zakharov (eds.), Springer- Verlag, Berlin 407-417 1993 



R-H APPROACH FOR CORREL. OF CHARACTERISTIC POLYNOMIALS 



35 



Johansson, K.: On fluctuations of eigenvalues of random herraitian matrices Duke Math. J. 91(1), 151-204 
1998 



Kanzieper, E ,: Replica field theories, Painleve transcendents and exact correlation functions, arxivxond 
/020724e| 



mat 



Keating, J., P., and Snaith, N., C: Random matrix theory and C(l/2 + it). Commun. Math. Phys. 214 , 
57-89 2000 

M.L. Mehta, Random Matrices, Academic Press, New York 1991 

Mehta, M., L., and Normand, J-M: Moments of the characteristic polynomial in the three ensembles of 
random matrices. J. Phys A: Math.Gen. 34, 4627-4639 2001 

Mirlin, A.D., and Fyodorov, Y.V.: Universality of level correlation function of sparse random matrices 
J.Phys.A: Math.Gen. 24 (10) 2273-2286 1991 

Nonnenmacher, S., and Zirnbauer, M., R.: Det-Det Correlations for quantum maps: Dual pair and saddle- 
point analyses J. Math. Phys. 43 (5), 2214-2240 2002 

Sagan, B.: The Symmetry Group Representations, Combinatorial Algorithms, and Symmetric Functions. 
Springer, New York 2000 



Splittorff, K and Verbaarschot, J., J., M.: 2002 Replica Limit of the Toda Lattice Equation arxivxond 



mat/0209594 



Szabo, R., J.: Microscopic spectrum of the QCD Dirac operator in three dimensions. Nucl Phys B 598 (1-2), 
309-347 2001 

Szego, G.: Orthogonal polynomials. American Mathematical Society, Colouium Publications, 23., Providence 
1975 

Tracy, C, A., and Widom, H.: Introduction to random matrices. Geometric and quantum aspects of intc- 
grable systems. Lecture Notes in Phys., vol. 424, Springer, Berlin 103-130 1993 

Tracy, C, A., and Widom, H.: Level spacing distribution and the Airy kernel. Comm. Math. Phys. 159, 
151-174 1994 

Tracy, C, A., and Widom, H.: Level spacing distribution and the Bessel kernel. Comm. Math. Phys. 161, 
289-309 1994 

Tracy, C, A., and Widom, H.: Fredholm determinants, differential equations and matrix models. Comm. 
Math. Phys. 163, 33-72 1994 

Zirnbauer, M., R.: Dual pairs in random matrix theory Talk given at LMS Workshop: Zeta Functions, 
Random Matrices and Quantum Chaos. September 13-14 2001 

Zirnbauer, M., R.: Random Matrices, Symmetry Classes, and Dual Pairs. Talk given at James H. Simons 
Workshop on Random Matrix Theory. Stony Brook, February 20 - 23 2002 



